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Abstract 

We compute the fourth virial coefficient of a binary nonadditive hard-sphere mixture over a wide 
range of deviations from diameter additivity and size ratios. Hinging on this knowledge, we build 
up a y expansion [B. Barboy and W. N. Gelbart, J. Chem. Phys. 71, 3053 (1979)] in order to trace 
the fluid-fluid coexistence lines which we then compare with the available Gibbs-ensemble Monte 
Carlo data and with the estimates obtained through two refined integral-equation theories of the 
fluid state. We find that in a regime of moderately negative nonadditivity and largely asymmetric 
diameters, relevant to the modelling of sterically and electrostatically stabilized colloidal mixtures, 
the fluid-fluid critical point is unstable with respect to crystallization. 
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INTRODUCTION 



A binary mixture of hard spheres is an established prototype model for understanding 
excluded-volume phenomena relevant to the equilibrium statistical mechanics of some com- 
plex fluids [l| • The model is characterized by the impenetrable diameters of the two species 
o\ and o"2 and by a cross diameter <j\i = (1/2) (o"i + (T 2 )(l + A), where the dimensionless 
quantity A accounts for deviations of the inter-species repulsive interactions from additiv- 
ity As is well known, in a real system the repulsive effects are associated with the overlap 
between the electron clouds of colliding molecules. However, even in a rare-gas mixture the 
Lorentz rule (A = 0) is not systematically satisfied [2]. Deviations from additivity are 
expected to become more and more pronounced with increasing pressures. In this respect, 
notwithstanding its crudeness, the hard-sphere model contains some basic features which can 
be used to sketch, for example, the phase behavior of binary fluids at very high pressures 
(~ 1 GPa) [3]. Deviations from additivity are also manifest in the "effective" interactions 
which characterize some complex polydisperse fluids such as colloid-polymer mixtures 4J. 

In general, a bidisperse mixture of hard spheres may exist in a surprisingly rich variety 
of both ordered and disordered phases whose formation is driven by entropic effects which 
depend on the diameter ratio q = cr 2 /cri (we conventionally identify the larger particles as 
species 1 in such a way that q < 1), on the deviation from additivity A, and on the mole 
fractions x% and x 2 = 1 — X\ of the two species, where X\ = iVi/(Aq + iV 2 ) [5|. For largely 
asymmetric size ratios (q < 0.3) and in the absence of nonadditive effects (A = 0), a mixture 
tends to phase separate, as a result of osmotic depletion effects, when the larger particles are 
sufficiently dilute {(J. However, Monte Carlo (MC) simulations show that in these conditions 
the demixing transition is preempted by the freezing of the larger particles 

The thermodynamic and phase-stability properties of a system turn out to be consid- 
erably affected by an even modest degree of nonadditivity. In a nonadditive hard-sphere 
(NAHS) mixture, the dominant mechanism underlying phase separation is different than in 
the additive case. First of all, the positive or negative nature of the nonadditivity is critical 
as to the type of thermodynamic behavior that is exhibited by the mixture. For A < the 
model shows a trend towards hetero-coordination and is able to reproduce the large com- 



positional fluctuations that are observed in some liquid and amorphous mixtures [8 



On the contrary, for A > even an equal-sized NAHS mixture exhibits phase separation 
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for high enough densities. This phenomenon is due to the extra repulsion between unlike 
spheres, which promotes the homo-coordination. The symmetric (equal-diameter) case has 



been investigated with a variety of theoretical methods 



HQ, 
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number of phenomenological equations of state have been proposed as well 
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22|. Fur- 



thermore, computer simulation studies have been performed for both negative and positive 



values of the nonadditivity parameter 
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25|, |26j. Correspondingly, the critical density 



as wel 



n 



as the universality class which the NAHS model belongs to have been determined 



So far, the asymmetric case has been studied much less than the symmetric one. The 
first simulation experiment of a NAHS mixture with similar diameters of the two species was 
carried out by Rovere and Pastore who traced the phase diagram by means of the Gibbs- 
ensemble Monte Carlo (GEMC) method 29|. Hamad performed Molecular Dynamics (MD) 
simulations and calculated the compressibility factor of a NAHS mixture for a number of 
values of q and A 30|. GEMC simulations of a highly asymmetric (q = 0.1) binary mixture 
show that fluid-fluid phase separation may occur on account of a rather small amount of 
nonadditivity 31]. as had been already conjectured by Biben and Hansen j^ . A reliable 
theoretical approach to study the phase diagram of NAHS mixture in a strongly asymmetric 
regime is based on the calculation of the effective potential between the larger particles. Such 
an effective potentia 
the smaller spheres 



is obtained upon integrating out the degrees of freedom associated with 



33 



35 



361 ] . Recently, we have shown, by means of integral-equation 



theories, that a stable fluid-fluid phase separation may occur for moderate size ratios and 
for a small amount of nonadditivity [sjjj. Another useful theoretical approach is based 
on the calculation of the residual multiparticle entropy complemented by MC numerical 
simulations [38]. This approach has been applied by Saija and Giaquinta to study the 
weakly asymmetric case (q > 0.75). 

Some efforts in the direction of studying NAHS mixtures by means of analytical theories 
have been also made by developing more refined equations of state (EOS) as compared with 
those obtained via a crude virial expansion. Schaink introduced an EOS whose derivation 
is reminiscent of the scaled particle theory and is valid in the symmetric as well as in the 
asymmetric case 39j. More recently, Santos and coworkers proposed an EOS that relies on 



the exact second and third virial coefficients and requires as an input the compressibility 
factor of the one-component system 401 ] . 
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In this paper we evaluate the fourth virial coefficient of a bidisperse NAHS mixture for 
a number of diameter ratios which span a wide range of accessible values and for —0.05 < 
A < 0.5. As far as we know, the only published virial coefficients beyond the third one are 
due to Vlasov and Masters [411 ] . However, these authors computed the virial coefficients of 
the mixture up to the sixth order for the size ratio q = 0.1 only and for just one positive 
value of the nonadditivity (A = 0.1). 

Even a reduced number of virial coefficients can be used to construct approximate expan- 
sions which turn out to be quite reliable over a wider density range than the original virial 
expansion [l| . It is in this conceptual framework that we aim at probing in this paper the re- 



liability of the so-called y expansion, originally proposed by Barboy and Gelbart (BG) [42 ]. 
We use this approach to calculate the fluid-fluid phase separation lines and compare the 
results with the available Gibbs-ensemble Monte Carlo data and with the results of some 
reliable integral-equation closures of the Ornstein-Zernike equations [371 ] . 

The paper is organized as follows: In Section II we present the numerical procedure used 
to calculate the fourth virial coefficient and the explicit expression of the BG EOS for a 
binary hard-sphere mixture. In Section III we present the phase diagram calculated for a 
number of size ratios. Some concluding remarks are finally given in Section IV. 

THE VIRIAL COEFFICIENTS AND THE BARBOY-GELBART EQUATION OF 
STATE 

The virial expansion can be written as: 

(3P = p + Bp 2 + Cp 3 + Dp 4 + ... (1) 

where Pis the pressure, (3 is the inverse temperature in units of the Boltzmann constant, and 
p = pi + p2, Pi {i = 1, 2) being the particle number density of the i-th species. In a mixture, 
at variance with the one-component case, the virial coefficients B,C, D . . . also depend on 
the relative concentration of the two species. In particular, the fourth-order coefficient reads: 

D = D llxl x{ + AD lll2 x\x 2 + §D ll22 x\x\ + AD l222 x x x\ + D 2222 x\ (2) 

In Eq.(T5]) the coefficients Dim and D 2222 can be calculated through the corresponding 
expression for a monodisperse fluid of particles with diameter <j\ or a 2 , respectively. The 
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cluster integrals Dm 2 and D n2 2 can be represented with the following four-point color 
graphs: 
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The open and solid circles identify in each graph particles belonging to species 1 and 2, 
respectively. Each bond contributes a factor to the integrand in the form of a Mayer step 
function. Space integration is carried out over all the vertices of the graph. We refer the 



reader for more details on the algorithm to 



43] 



The coefficient -D1222 is obtained from the general expression of D1112 by merely inter- 
changing the larger with the smaller particles and vice-versa. The Monte Carlo results for 
the composition-independent coefficients D a ^s (with a, /3,7, S = 1,2) are given in tables 
Ul lVIl over the whole range of q for A < 0.5. As q — >• 1, all the coefficients tend to the value 
2.6362... (in units of af) of the one-component fluid. In the opposite limit of q — > the 
partial coefficient D U12 approaches the exact result l/4(7rcrf/6) [43j. The knowledge of the 
fourth virial coefficient allows one to derive the BG expansion up to the fourth order upon 
expanding the EOS in powers of the variables z/i = Pi/(1 — r]), where r\ = tc p x a\ / '6 + tt P2<j\ / '6: 
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are determined so as to ensure the low-density expansion of Eg. (151) to coincide with the 
fourth-order virial expansion reported in Eq. ([!]). The quantities B a p and C al 3-y are the 



Tl ' 
111- 



In the 



second and third virial coefficients, respectively, which are available analytically 
limit A — > 0, the first three terms of Eq.(jSD, which correspond to the third order expansion 
of the BG EOS, reduce to the compressibility equation-of-state of the Percus-Yevick theory. 
Upon rearranging the above relations, we can write: 

(i - >?)' 

where 

Ai(x) = 1 

A 2 {x) = x\Au + 2xi(l - xx)A 12 + (1 - x 1 ) 2 A 22 
A 3 (x) = x\Am + 3a;i(l - x x )A 112 + 3x(l - x x ) 2 A 122 + (1 - x x f A 222 
A 4 (x) = xt A U11 +Axl(l-x 1 )A 1U2 + 6x1(1 -x 1 ) 2 A U22 + Ax 1 (l-x 1 ) 3 A 1222 
+ (1 - xi) 4 A 2222 

As we shall see in the following Section, the inclusion of the fourth virial coefficient 
definitely improves the performance of the BG EOS over the third-order one. A more 
convincing test of the reliability of the BG EOS including the fourth virial coefficient is 
provided by the comparison of the predictions for the fluid-fluid phase separation threshold 
with the numerical simulation data. In order to trace the phase-coexistence curve, we need 
the Gibbs free energy. We start by calculating the Helmholtz free energy: 



PA r ( z -w 



N Jo V 



(PA 3 ) " 1] + X> H*k) - ln(l (10) 



fe=i fe=i 



where Z = f3P/ p is the compressibility factor; the Gibbs free energy then follows as 44. |45| : 

^ = ^ + Z. (11) 

The fluid-fluid coexistence curves are obtained through the common-tangent construction 
of the Gibbs free energy at constant pressure. 



RESULTS AND DISCUSSION 



We present in Figs. [T] and [2] the compressibility factor calculated for different values of 
the size ratio and of the nonadditivity parameter. The fourth-order BG expansion (BG4) of 



the EOS is in better agreement with the simulation data provided by Hamad [30J than the 
corresponding third-order expansion (BG3), expecially for large packing fractions. However, 
as shown in Fig. [21 this improvement is not systematic. 

Species with not too dissimilar sizes 

We compare in Fig. [3] our results for the phase coexistence thresholds with those obtained 
by Rovere and Pastore for q = 0.8333 and A = 0.182 in their GEMC simulations j^. Un- 
fortunately, we could not report in the figure the simulation critical parameters explicitly 
because the coexistence curve becomes rather flat around there and a reasonable determi- 
nation of them would typically require very long computations [29]. As shown in the top 
panel of Fig. [3J, both BG estimates of the EOS appear to underestimate the critical pressure, 
which is located at the minimum of the phase coexistence curve, by about 20%. However, on 
account of the large fluctuations that affect the system on approaching the critical point, it 
may be possible that the pressure measured inside the two simulation boxes is not sufficiently 
stable against the value imposed from outside. This fact might be partially responsible for 
the discrepancy observed between the GEMC data and the BG estimates. The growth of 
fluctuations near the critical point is apparent in the bottom panel of Fig. [3] where we report 
the error bars associated with the GEMC results for the coexistence densities plotted as a 
function of the relative composition. In fact, in the density-concentration plane, the present 
theoretical predictions fall inside the simulation error bars around the critical point. Overall, 
the BG4 critical pressure and density show a better agreement with the GEMC simulation 
results as compared with their BG3 counterparts. 

We also show the behavior of the analytical EOS for size ratios q = 0.75 (Fig. H]) and 
q = 0.6 (Fig. [5]) and for two values of the nonadditivity parameter A = 0.05,0.1 which 



are particularly meaningful for Helium-Xenon mixtures at high pressure [37|]. In this case 
the critical point falls at medium-high packing fractions and we found no simulation data 
in the literature to compare with. Hence, we compare the BG3 and BG4 estimates of 
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the phase separation curves with the predictions of two integral-equation theories 46JI, i.e., 
the modified hypernetted-chain (MHNC) [47J and the Rogers- Young (RY) [48] closures of 
the Ornstein-Zernike equations. In these approximate theories two internal parameters are 
adjusted so as to ensure the equality of the two osmotic compressibilities that are respectively 
obtained via number fluctuations or upon differentiating the virial pressure 37]]. We show 
in Fig. H] (upper panel) the binodal line for q = 0.75 and A = 0.05. The current numerical 
estimates for the critical mole fraction and for the reduced pressure are reported in Table [VTTl 
It is evident that the BG4 critical pressure is in better agreement with the corresponding 
integral-equation estimates. The differences between the BG3 and BG4 expansions are 
smaller for A = 0.1 (see Fig.HJ bottom panel). Even in this case, the BG4 critical parameters 
are in better agreement with the corresponding MHNC and RY estimates. As for the 
shape of these curves, we note that the values of the constant-pressure Gibbs free energy 
obtained from integral-equation theories was fitted through a fourth-order polynomial curve 
(see Fig. El upper panel) [371. Overall, the resulting agreement with the current results for the 
Gibbs free energy is good. However, the implementation of the Maxwell construction requires 
the concentration derivative of the Gibbs free energy (see Fig. [6]). The discrepancies observed 
in the shapes of the phase-separation curves obtained with the two methods discussed above 
(EOS expansions and integral-equation theories) may be likely due to the estimate of this 
latter quantity. In fact, we suspect that the derivative of the polynomial curve used to fit 
the integral-equation data for the Gibbs free energy probably is not flexible enough as one 
is lead to argue after observing the shapes of the corresponding BG3 and BG4 curves. 

The BG4 fluid-fluid critical point is in better agreement with the corresponding integral- 
equation estimate also for q = 0.6 (see Fig. [5]). In fact, for A = 0.05 the BG4 critical 
pressure, that is reported in Table IVHH falls in between the MHNC and RY values as for 
q = 0.6 (see Table |VH| . As far as the critical concentration is concerned, in passing from 
q = 0.75 to q — 0.6, the BG4 estimate moves to lower x\ values in a more distinct way than 
its MHNC and RY courterparts. Finally, for larger values of the nonadditivity parameter 
(A = 0.1), the differences in the critical parameters estimated with these approaches smooth 
over, as previously shown for q = 0.75 (see Table |YTH| . 
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Highly asymmetric regime 



We start commenting the results obtained for q = 0.1 and nonadditivity A = 
0.2,0.3,0.4,0.5. We compare in Fig. [7] the fluid-fluid coexistence curves obtained by means 
of the BG3 and BG4 estiamtes of the EOS with the GEMC simulations reported in [31]. As 
already observed for larger size ratios, the BG4 curve shows a better agreement with the 
simulation data in the region around the critical point, and this improvement is more evident 
for low values of the nonadditivity parameter. On the other side, for larger values of A, the 
BG3 and BG4 estimates become very close to each other, as was also the case of mixtures 
with not too dissimilar species. This effect is probably related to the comparatively lower 
packing fractions of the demixing region that is observed for larger nonadditivities. In fact, 
the mixture is less correlated at low densities and one may then expect that the differences 
between the fourth and the third-order estimates of the EOS tend to disappear. Figure [8] 
shows the critical pressure and packing fraction plotted as a function of the nonadditivity 
parameter. We also observed a crossover in the BG4 critical pressure versus its BG3 value 
upon further decreasing A: viz., the BG4 values get lower than the BG3 ones. This obser- 
vation is particularly significant because the BG3 critical pressure tends to diverge in the 
limit of additive mixtures (A = 0), whilst the BG4 value apparently remains finite. 

As pointed out by many authors, the regime of low or even negative values of the nonad- 
ditivity parameter is interesting for the application of the model to sterically and electrostat- 
ically stabilized colloidal mixtures. In particular, Louis and coworkers used the Alexander- 
de Gennes theory with the Derjaguin approximation in order to estimate the value of the 
nonadditivity parameter for a polymethylmetrylacrilate mixture stabilized by a poly- 12- 
hydroxysteric acid brush {33]. They obtained A = —0.01 for q = 0.1075. We present in the 
bottom panel of Fig. [9] the BG4 fluid-fluid binodal line with its critical point for a mixture 
with size ratio q = 0.1 and compare it with the fluid-solid binodal line as reported in 33]. 
For negative values of A, the demixing transition should be preempted by the formation of 



the solic 
only 



36 



Dhase because a stable demixing is expected to occur in very asymmetric mixtures 



49j. In ref. [33J, the adopted crystal structure in order to evaluate the fluid-solid 



binodal, is an FCC lattice of large spheres (permeated by a fluid of small spheres). Actually, 
we find that demixing occurs for very high values of the larger-particles packing fraction 
(see the inset of Fig. [9]) and for very low values of the smaller-particles packing fraction. 
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The corresponding BG4 critical point, with coordinates rji = 0.67 and 7/2 = 0.0014, falls well 
inside the region of stability of the solid phase. Louis and coworkers also estimated that 
for positive values of the nonadditivity parameter the critical point remains metastable up 
to A « 0.2 [33|. In agreement with their prediction, for A = 0.25 the fluid-fluid critical 
point, located by both the BG3 and BG4 expansions, falls below the freezing line (see the 
top panel of Fig. [9]) . 



CONCLUDING REMARKS 



In this paper we discussed the fourth virial coefficient of binary hard-sphere mixtures, 
which has been computed for different values of the nonadditivity parameter and of the size 
ratio. We have shown that its use within the y expansion originally proposed by Barboy 
and Gelbart improves over the description of the fluid-fluid phase separation provided by the 
same expansion truncated at the third order, the more so for small values of the size ratio 
and of the nonadditivity parameter. In particular, we found that the pressure calculated at 
the critical point through the fourth-order expansion does not diverge in the limit of small 
values of the nonadditivity parameter and of largely asymmetric size ratios. This feature 
allowed us to predict that, within this approximation, the critical point is unstable with 
respect to crystallization in the regime relevant to colloidal mixtures. 
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TABLE I: Fourth-order partial virial coefficients as a function of the size ratio for A = 0.05. The 
numerical values are given in units of a\. The error on the last significant figure is enclosed in 
parentheses. 



Q 


-Dim 


-C1122 


-C1222 


0.05 


0.0835(3) 


5.53(1)-10~ 5 


2.064(1)-10~ 8 


0.10 


0.1201(4) 


5.271(9)-10- 4 


1.3981(6)-10- 6 


0.15 


0.1661(5) 


2.073(4)-10- 3 


1.6840(8)-10- 5 


0.20 


0.2225(6) 


5.67(1)-10~ 3 


9.988(6)-10~ 5 


0.25 


0.2900(7) 


1.267(2)-10~ 2 


4.019(2)-10~ 4 


0.30 


0.3697(8) 


2.485(4)-10~ 2 


1.2645(8)-10- 3 


0.40 


0.569(1) 


7.46(1)-10- 2 


7.857(6)-10- 3 


0.50 


0.828(1) 


0.1806(4) 


3.303(3)-10- 2 


0.60 


1.153(2) 


0.3807(8) 


0.1084(1) 


0.70 


1.554(2) 


0.728(1) 


0.2990(3) 


0.75 


1.784(2) 


0.977(1) 


0.4726(3) 


0.80 


2.036(3) 


1.292(2) 


0.7269(8) 


0.85 


2.309(2) 


1.683(2) 


1.0914(9) 


0.90 


2.606(3) 


2.166(4) 


1.604(2) 


0.95 


2.927(4) 


2.755(6) 


2.311(3) 
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TABLE II: Fourth-order partial virial coefficients as a function of the size ratio for A = 0.1 (see 
also the caption of Table I) . 



q 


-Dim 


-C1122 


-C1222 


0.05 


0.1201(4) 


5.37(1)-10~ 5 


2.3824(8)-10~ 8 


0.10 


0.1686(7) 


6.15(2)-10- 4 


1.620(1)-10- 6 


0.15 


0.2284(6) 


2.505(4)-10~ 3 


1.959(1)-10~ 5 


0.20 


0.3007(7) 


6.94(1)-10~ 3 


1.1664(6)-10~ 4 


0.25 


0.3863(8) 


1.560(2)-10~ 2 


4.710(3)-10~ 4 


0.30 


0.4865(9 


3.073(6)-10~ 2 


1.486(1)-10- 3 


0.40 


0.735(1) 


9.26(2)-10- 2 


9.298(7)-10~ 3 


0.50 


1.053(2) 


0.2249(4) 


3.934(3)-10- 2 


0.60 


1.451(2) 


0.4750(8) 


0.1298(1) 


0.70 


1.938(3) 


0.909(2) 


0.3602(3) 


0.75 


2.216(2) 


1.220(2) 


0.5710(4) 


0.80 


2.519(3) 


1.613(3) 


0.8810(9) 


0.85 


2.849(2) 


2.103(2) 


1.326(1) 


0.90 


3.206(4) 


2.704(5) 


1.953(2) 


0.95 


3.591(4) 


3.442(6) 


2.822(3) 
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TABLE III: Fourth-order partial virial coefficients as a function of the size ratio for A = 0.2 (see 
also the caption of Table I) . 



Q 


-D1112 


-C1122 


-D1222 


0.05 


0.2226(6) 


-2.277(2)-10~ 4 


3.116(1)-10- 8 


0.10 


0.3007(7) 


1.15(2)-10- 4 


2.134(1)-10~ 6 


0.15 


0.3945(8) 


2.022(6)-10~ 3 


2.595(1)-10~ 5 


0.20 


0.506(1) 


7.10(2)-10~ 3 


1.5554(8)-10~ 4 


0.25 


0.637(1) 


1.758(4)-10~ 2 


6.320(4)-10~ 4 


0.30 


0.787(1) 


3.647(7)-10- 2 


2.006(1)-10- 3 


0.40 


1.154(2) 


0.1150(2) 


1.2694(9)-10 -2 


0.50 


1.617(2) 


0.2879(5) 


5.429(4)-10- 2 


0.60 


2.190(3) 


0.616(1) 


0.1809(2) 


0.70 


2.881(4) 


1.188(2) 


0.5068(5) 


0.75 


3.274(3) 


1.600(2) 


0.8071(5) 


0.80 


3.702(4) 


2.118(4) 


1.250(1) 


0.85 


4.166(3) 


2.765(4) 


1.890(1) 


0.90 


4.664(5) 


3.562(6) 


2.795(3) 


0.95 


5.204(6) 


4.528(8) 


4.056(4) 
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TABLE IV: Fourth-order partial virial coefficients as a function of the size ratio for A = 0.3 (see 
also the caption of Table I) . 



Q 


-Dim 


-C1122 


-C1222 


0.05 


0.3697(8) 


-2.0687(2)-10~ 4 


3.984(1)-10~ 8 


0.10 


0.4865(9) 


-3.736(3)-10- 3 


2.744(1)-10~ 6 


0.15 


0.625(1) 


-4.857(7)-10~ 3 


3.357(1)-10~ 5 


0.20 


0.787(1) 


-3.90(2)-10~ 3 


2.022(1)-10~ 4 


0.25 


0.974(2) 


1.21(5)-10~ 3 


8.261(5)-10" 4 


0.30 


1.188(2) 


1.339(9)-10- 2 


2.635(2)-10~ 3 


0.40 


1.705(2) 


7.46(3)-10~ 2 


1.683(1)-10 -2 


0.50 


2.350(3) 


0.0.2219(7) 


7.259(5)-10~ 2 


0.60 


3.139(4) 


0.514(2) 


0.2439(2) 


0.70 


4.085(5) 


1.034(3) 


0.6886(6) 


0.75 


4.622(4) 


1.410(2) 


1.1001(7) 


0.80 


5.203(6) 


1.188(5) 


1.711(2) 


0.85 


5.833(4) 


2.477(4) 


2.594(2) 


0.90 


6.506(7) 


3.204(9) 


3.849(4) 


0.95 


7.233(7) 


4.09(1) 


5.602(5) 
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TABLE V: Fourth-order partial virial coefficients as a function of the size ratio for A = 0.4 (see 
also the caption of Table I) . 



Q 


-D1112 


-D1122 


-D1222 


0.05 


0.569(1) 


-9.0438(3)-10~ 3 


5.004(1)-10~ 8 


0.10 


0.734(1) 


-1.7721(2)-10~ 2 


3.462(2)-10~ 6 


0.15 


0.929(1) 


-3.002(7)-10~ 2 


4.256(4)-10~ 5 


0.20 


1.153(2) 


-4.597(3)-10~ 2 


2.575(1)-10~ 4 


0.25 


1.412(2) 


-6.544(6)-10" 2 


1.0564(6)-10- 3 


0.30 


1.705(2) 


-8.75(1)-10- 2 


3.383(2)-10~ 3 


0.40 


2.406(3) 


-0.1377(4) 


2.178(1)-10~ 2 


0.50 


3.274(4) 


-0.1898(9) 


9.461(7)-10~ 2 


0.60 


4.327(5) 


-0.236(2) 


0.3199(2) 


0.70 


5.584(6) 


-0.280(4) 


0.9088(7) 


0.75 


6.295(5) 


-0.299(3) 


1.4564(9) 


0.80 


7.062(7) 


-0.314(6) 


2.270(2) 


0.85 


7.890(6) 


-0.346(6) 


3.452(2) 


0.90 


8.782(9) 


-0.39(1) 


5.137(5) 


0.95 


9.729(9) 


-0.46(1) 


7.498(7) 
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TABLE VI: Fourth-order partial virial coefficients as a function of the size ratio for A = 0.5. See 
also caption of table I. 



q 


-Dim 


-D1122 


-C1222 


0.05 


0.828(1) 


-2.91908(3)-10~ 2 


6.180(1)-10~ 8 


0.10 


1.053(2) 


-5.6318(3)-10" 3 


4.295(2)-10~ 6 


0.15 


1.316(2) 


-9.806(1)-10- 2 


5.303(4)-10- 5 


0.20 


1.617(2) 


-0.15883(3) 


3.220(2)-10~ 4 


0.25 


1.962(3) 


-0.24429(7) 


1.3258(7)-10~ 3 


0.30 


2.350(3) 


-0.3605(1) 


4.263(2)-10~ 3 


0.40 


3.275(4) 


-0.7246(5) 


2.761(2)-10~ 2 


0.50 


4.410(5) 


-1.346(1) 


0.12067(8) 


0.60 


5.781(6) 


-2.378(2) 


0.4103(3) 


0.70 


7.409(7) 


-4.064(4) 


1.1717(9) 


0.75 


8.326(6) 


-5.245(4) 


1.882(1) 


0.80 


10.314(9) 


-6.718(7) 


2.941(2) 


0.85 


10.379(7) 


-8.578(7) 


4.482(3) 


0.90 


11.52(1) 


-10.88(1) 


6.686(6) 


0.95 


12.74(1) 


-13.73(2) 


9.777(8) 
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TABLE VII: Critical molar fractions and reduced pressures for q = 0.75. 



A 




Xi 


{5Po\ 


0.05 


BG3 


0.34 


12.21 




BG4 


0.36 


10.34 




MHNC 


0.44 


9.72 




RY 


0.43 


10.5 


0.1 


BG3 


0.36 


4.73 




BG4 


0.37 


4.54 




MHNC 


0.455 


4.1 




RY 


0.45 


4.12 



TABLE VIII: Critical molar fractions and reduced pressures for q = 0.60. 



A 






f3Paf 


0.05 


BG3 


0.241 


16.85 




BG4 


0.255 


14.15 




MHNC 


0.405 


13.7 




RY 


0.36 


15.05 


0.1 


BG3 


0.27 


6.55 




BG4 


0.28 


6.30 




MHNC 


0.41 


5.74 




RY 


0.41 


6.25 
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FIG. 1: Compressibility factor plotted as a function of the total packing fraction for q = 1/3, 
x% = 0.5: left panel, A = —0.05; right panel, A = 0.05. The continuous and dashed lines represent 
the BG4 and BG3 estimates, respectively. The solid circles are the MD results by 30J]. The cor- 



responding fourth-order partial virial coefficients calculated for A = —0.05 are: D1112=0. 2304(6), 
D1122=2.037(5) • 10~ 2 , D1222=l. 719(1) • 10~ 3 . 
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FIG. 2: Compressibility factor plotted as a function of the mole fraction for A = 0.2: left panel, 



q = 0.25; right panel, q = 0.5. Stars: BG4; triangles: BG3; circles: MD results 



30] 
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FIG. 3: Fluid-fluid coexistence in the pressure-composition plane (top panel) and density- 
composition plane (bottom panel) plotted for q = 0.8333 and A = 0.182. The continuous and 
dashed curves represent the BG4 and BG3 binodal lines, respectively. The critical point is identi- 
fied with solid (BG4) and open (BG3) diamonds. Open circles, with the corresponding error bars 
in the density-composition plane, represent the GEMC simulation data [29(]. The corresponding 
fourth-order partial virial coefficients are: D1112=3. 752(5), D1122=2.457(3), D1222=1.459(l). 
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FIG. 4: Fluid-fluid coexistence in the pressure-composition plane for q = 0.75: top panel, A = 0.05; 
bottom panel, A = 0.1. The continuous and dashed curves represent the BG4 and BG3 binodal 
lines, respectively. The critical point is identified with solid (BG4) and open (BG3) diamonds. The 



dotted and dot-dashed lines represent the RY and MHNC predictions, respectively 



3J. 
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FIG. 5: Total Gibbs free energy (top panel) and its first derivative (bottom panel) plotted as a 
function of the mole fraction for q = 0.75, A = 0.05 and f3Paf = 11.38. Dashed line: BG3; 
continuous line: BG4; solid circles: RY closure; dot-dashed line: fourth-order polynomial fit of the 
RY data. 



24 




FIG. 6: Fluid-fluid coexistence in the pressure-composition plane for q = 0.6: top panel, A = 0.05; 
bottom panel, A = 0.1. The continuous and dashed curves represent the BG4 and BG3 binodal 
lines, respectively. The critical point is identified with solid (BG4) and open (BG3) diamonds. The 



dotted and dot-dashed lines represent the RY and MHNC predictions, respectively 



3J.. 
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FIG. 7: Fluid-fluid coexistence in the pressure-composition plane for q = 0.1. The continuous and 
dashed curves are the BG4 and BG3 binodal lines, respectively. The critical point is identified 
with solid (BG4) and open (BG3) diamonds. The open symbols represent the GEMC simulation 
data [3l| obtained for A = 0.5 (circles), 0.4 (squares), 0.3 (diamonds), 0.2 (triangles). 
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FIG. 8: Reduced pressure (top panel) and total packing fraction (bottom panel), calculated at the 
critical point for a nonadditive hard-sphere mixture with size ratio q = 0.1, plotted as a function 
of the nonadditivity parameter A: open diamonds, BG3; solid diamonds, BG4; stars, simulation 
data [311 ]. 
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FIG. 9: Phase diagram calculated for q = 0.2 and A = 0.25 (top panel) and for q = 0.1075 and 
A = —0.01 (bottom panel). Dot-dashed curves: fluid-solid binodal line from [331 ] : the continuous 
and dashed curves represent the BG4 and BG3 binodal lines, respectively. The critical point is 
identified with solid (BG4) and open (BG3) diamonds. The continuous lines plotted in the inset 
represents the BG4 fluid-fluid binodal line. The orresponding fourth-order partial virial coefficients 
are: D1112=7.98(l) • 10~ 2 , D1122=4.840(4) • 10" 4 , D1222=1.8029(7) • 10~ 6 . 
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